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ABSTRACT 

Aspects  of  crack  growth  in  an  elastic-plastic  material  under  quasi-static, 
steady  conditions  are  investigated  for  plane  problems.  A  path-independent  line 
integral  is  identified  for  the  steady  problem  which  generalizes  the  J-integral  of 
deformation  theory  plasticity  to  incremental  theories  of  plasticity.  Implications 
of  the  integral  for  small  scale  yielding  crack  growth  are  discussed.  A  model 
study  is  made  of  the  effect  on  the  growth  process  of  the  residual  plastic  wake 
left  behind  the  advancing  crack.  The  influence  of  the  wake  on  the  stabilization 
of  crack  growth  is  discussed. 


1.  INTRODUCTION 

Irreversibility  effects  which  tend  to  stabilize  crack  growth  arise  from  two 
related  sources:  nonproportional  plastic  deformation  in  the  active  plastic  zone 
and  elastic  unloading  which  leaves  a  wake  of  residual  plastic  deformation  behind 
the  advancing  crack  tip.  McClintock  [1,  2]  studied  quasi-static  crack  growth 
under  anti-plane  shear  in  an  eiastic-perfectly  plastic  material.  Of  particular 
significance  is  his  discovery  that  the  strains  in  the  active  plastic  zone  depend 
logarithmically  on  the  distance,  r  ,  from  the  crack  tip  as  the  tip  is  approached. 
The  corresponding  behavior  for  a  stationary  crack  in  a  perfectly  plastic  material 


This  work  was  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research 
under  Grant  AFOSR-73-2476,  in  part  by  the  Advanced  Research  Projects  Agency 
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is  far  more  singular  with  strain  varying  like  1/r  .  McClintock  was  able  to 
estimate  the  strain  field  in  the  region  ahead  of  the  crack  for  small  scale  yielding. 
Then  employing  a  fracture  criterion  based  on  the  attainment  of  a  critical  strain 
at  some  characteristic  ahead  of  the  crack,  he  showed  that  stable  crack  growth 
should  be  expected  for  any  reasonable  choice  of  critical  strain  and  characteristic 
distance.  Later,  Chitaley  and  McClintock  [3]  carried  out  a  detailed  numerical 
analysis  of  the  small  scale  yielding  problem  under  steady  growth  conditions  and 
verified  the  major  features  of  the  more  approximate  analysis  given  previously  in 
[1.  2]. 

Rice  [4,  5]  has  shown  that  a  logarithmic  dependence  of  the  strains  on  r  can 
also  be  expected  for  a  growing  crack  in  an  elastic-perfectly  plastic  material  for 
both  plane  stress  and  plane  strain.  More  complete  analyses  of  the  plane  problems, 
analogous  to  those  cited  above  for  anti-plane  shear,  have  not  been  achieved.  None¬ 
theless,  because  of  the  weak  logarithmic  singularity  in  the  strains,  it  seems 
probable  that  such  an  analysis  would  again  predict  stable  crack  growth  for  the 
plane  problems.  If  true,  this  would  be  cause  for  some  concern  since  under  conditions 
approximating  plane  strain  little  or  no  stable  crack  growth  is  frequently  observed 
in  small  scale  yielding  tests.  Since  the  stress  history  of  each  material  point 
swept  by  the  advancing  plastic  zone  is  distinctly  nonproportional  [3,  5],  it  is 
conceivable  that  the  smooth  yield  surface  of  the  plasticity  theory  used  in  the 
singularity  analysis  overly  restricts  the  plastic  flow  and  leads  to  a  weaker 
singularity  in  the  strains  than  might  otherwise  be  found.  For  similar  reasons, 
a  smooth  yield  surface  is  thought  to  be  inadequate  for  describing  plastic  flow  in 
the  bifurcation  analysis  of  plastic  buckling. 

In  this  paper  the  influence  of  the  residual  plastic  wake  on  the  extending  crack 
will  be  modeled  for  steady  conditions.  The  model  is  similar  in  spirit  to  the 
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Dugdale-Barenblatt  model  for  a  stationary  crack  in  a  thin  sheet.  The  model  is  not 
intended  to  be  complete  since  irreversibility  effect?,  associated  with  nonproportional 
loading,  such  as  those  discussed  above,  are  not  included.  Instead,  the  purpose  of 
the  model  is  to  complement  the  above  mentioned  studios  to  illustrate  that  the  wake 
by  itself  can  have  a  significant  effect  on  stabilizing  crack  propagation.  The  model 
is  formulated  for  plane  stress  but  some  difference  between  plane  stress  and  plane 
strain  can  be  inferred  with  respect  to  the  effect  of  the  wake  on  growth. 

Before  introducing  the  model,  a  general  derivation  of  an  energy  balance  relation 
is  made  for  steady  propagation  of  cracks  in  elastic-plastic  materials.  A  path- 
independent  integral  is  Identified  which  generalizes  the  J-integral  of  Rice  [4]  to 
hold  for  arbitrary  material  behavior.  Our  development  is  closely  related  to  the  work 
of  Cherepanov  [6,  7].  The  integral  is  used  to  show  the  relation  between  the  residual 
internal  energy  in  the  wake  and  the  elastic  stress  intensity  factor  for  small  scale 
yielding.  This  information  is  then  used  in  constructing  the  model. 


2.  A  PATH-INDEPENDENT  INTEGRAL  FOR  STEADY  CRACK  GROWTH  IN  THE  PLANE 
Let  the  coordinates  x^  and  X£  translate  with  the  crack  tip  and  be  centered 
as  in  Fig.  1.  Crack  growth  occurs  in  the  x^  direction;  let  (X^  *  a,  X£  *  0)  be 
the  position  of  the  crack,  tip  with  respect  to  some  fixed  coordinate  system  (X^^^)  . 
The  crack-tip  coordinate  a  is  assumed  to  increase  monotonically  and  can  be  taken 
as  the  time-like  variable  in  the  quasi-static  growth  process.  If  steady  conditions 
hold,  either  globally  or  in  some  vicinity  of  the  crack  tip,  then  the  stresses,  strains 
and  displacements  can  be  expressed  as  functions  of  the  translating  coordinates  and 

X2  independent  of  a  .  The  following  relations  among  the  derivatives  hold: 


(1) 


Small  strain  plasticity  is  considered  in  which  the  equilibrium  equations  and 
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-he  strain-uisplaceraent  relations  are  taken  to  be  linear.  Thermal  effects  and 
inertia  are  neglected.  Let  W  be  the  internal  energy  density  at  any  material 
point,  i.e.,  the  stress  work  density 

E. 


W 


ij 


aijd£ij 


(2) 


Here  the  notation  of  Rice  [4]  is  used;  however,  we  assume  only  that  the  behavior 
of  the  material  is  time-independent  and  we  do  not  limit  consideration  to  a 
deformation  theory  of  plasticity.  The  integral  in  (2)  will  depend  in  general  on 
the  stress  history  at  the  material  point  in  question. 

Let  T  be  any  closed  curve  in  the  plane  (not  enclosing  the  tip)  with  enclosed 
area  Ap  which  translates  with  the  crack  tip  without  deforming.  Let  $p  be  the 
internal  energy  per  unit  thickness  enclosed  within  T  so  that 


-Jw 


dA 


(3) 


Under  steady  conditions  <J>p  is  constant  and  thus 


0  - 


d#r 

da 


dA  + 


W  n^  ds 


(4) 


where  n^  is  the  outward  unit  normal  to  T  and  ds  is  the  length  element.  Next 


note  that 


l  (£]t  dA  ■  1  41  ■  -\ 


dA 


-I 


Vu’,1  “ '  ‘ 


°ij”jui,i  ds 


(5) 
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where  a..  -  0  and  e,,  -  ^(u  +  u  )  have  also  been  used,  and  where 

i  j  9  j  1 3  ^  1  >  J  J  » 1 

(  )  =  3(  )/3x  .  Combining  (4)  and  (5)  gives,  for  any  closed  contour, 

J  1  *• 

(W  n^T^  1>ds  “  0 


where  Tj 


Vi 


The  path-independence  of  this  integral  can  be  demonstrated  directly  by  applying 
Green's  theorem  to  the  first  term  in  (6), 


JW  n^  ds  -  j  W  ^  dA 

r  Ar 


(7) 


and  then  noting  that  (5)  applies.  The  derivation  assumes  there  are  no  singularities 
inside  T  .  But  it  still  holds  when  discontinuities  in  the  second  derivatives  of 
the  displacements  and  in  the  first  derivatives  of  the  stresses  occur  across 
contours  contained  within  T  ,  as  can  be  expected  across  the  instantaneous  boundary 
between  the  plastically  loading  and  elastically  unloaded  regions. 

Each  of  the  two  terms  in  (6)  has  a  simple  physical  interpretation.  The  first 
is  the  rate  per  unit  crack  extension  at  which  stress  work  passes  through  the  contour 
T  which  is  translating  with  the  tip.  The  second  is  the  rate  of  work  per  unit  crack 
extension  done  by  the  tractions  acting  on  T  on  the  material  within.  Since,  the 
total  stress  work  within  T  is  constant  in  the  steady  state,  the  sum  of  these  two 
terms  must  vanish. 

Now  consider  a  closed  contour  such  as  that  shown  in  Fig.  1  where 
T  -  ri+  r2+  .  The  contributions  along  the  crack  face,  and  ,  vanish 

identically.  Thus  the  integral  on  equals  the  integral  on  . 


That  is,  the 
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integral  on  any  contour  ,  such  as  or  ,  which  encloses  the  crack  tip 

is  the  same.  Let 


J 


s 


(Wni-Vlii)ds 


(8) 


c 

where  the  superscript  s  is  to  emphasize  that  it  is  restricted  to  steady  conditions. 


The  similarity  in  notation  between  (8)  and  the  integrals  discussed  by  Eshelby 
[8]  and  Rice  [4]  is  deceptive  since  (8)  is  path-independent  for  any  material  as 
long  as  steady  conditions  pertain.  In  general,  (8)  is  not  path-independent  when 
crack  propagation  is  not  steady.  In  contrast,  the  integral  discussed  by  Eshelby 
and  Rice  is  restricted  to  nonlinear  elastic  solids  (e.g.,  a  deformation  theory  of 
plasticity)  but  is  not  restricted  to  steady  conditions. 

We  now  argue  that  J8  =  0  for  steady  growth  of  a  line  crack  in  an  elastic- 

plastic  material  under  the  small  strain  assumptions  for  which  (8)  was  derived.  For 

elastic-perfectly  plastic  problems  involving  a  line  crack  in  anti-plane  shear,  plane 

strain  or  plane  stress  it  has  been  noted  that  the  singularity  in  the  strains  (and 

more  generally  the  displacement  gradients)  at  the  crack  tip  is  only  logarithmic. 

It  follows  that  the  singularity  of  the  integrand  in  (8)  is  also  logarithmic. 

Consequently,  if  F^  is  taken  as  a  small  circle  or  radius  r£  centered  at  the  tip 

so  that  ds  ■  rc  d0  ,  then  by  letting  r£  -*•  0  it  is  immediately  seen  that  J8  *  0  . 
s 

Physically,  J  *  0  implies  that  no  energy  feeds  out  at  the  advancing  tip  of  the 
line  crack  in  the  elastic-plastic  material  as  it  does  in  a  linear  or  nonlinear  elastic 
material.  This  same  procedure  has  quite  a  different  conclusion  when  the  material  is 
nonlinear  elastic.  Then,  it  is  known  that  J  +  0  and  the  above  procedure  implies 
that  the  integrand  in  (8)  must  have  a  1/r  singularity,  which  has  been  verified 
for  a  class  of  nonlinear  materials  in  [9,  10]. 
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No  information  is  available  on  the  singular  behavior  of  the  stress  and  strain 
at  the  crack  tip  in  a  strain  hardening  material  for  nonstationary  cracks.  However 
it  would  not  be  expected  that  a  small  degree  of  strain  hardening  could  lead  to  a 
dramatically  stronger  singularity  in  W  than  that  occurring  in  the  perfectly  plastic 
case.  In  fact,  we  suspect  that  any  amount  of  irreversible  plastic  flow  in  the  steady 
problem  will  lead  to  a  singularity  in  W  which  is  weaker  than  1/r  so  that  Js  *  0  . 

The  path-independent  integral  (8)  is  closely  related  to  an  integral  derived  by 
Cherepanov  [6],  although  in  several  respects  our  treatment  differs  from  his.  Most 

g 

Importantly,  our  conclusion  that  J  vanishes  for  steady  growth  is  in  direct 
contradiction  to  his  conclusion.  As  in  this  paper,  Cherepanov  considers  general 
material  behavior.  However,  he  does  not  restrict  consideration  to  steady  crack  growth. 
Following  a  series  of  physical  and  mathematical  arguments,  he  arrives  at  the 
propagation  criterion  (in  the  present  notation  and  neglecting  thermal  and  inertial 
effects  which  are  also  considered  in  [6}): 

(W  n1-T1ui>1)ds  -  2y  (9) 

c 


Consider  the  following  argument.  Suppose  irreversible  plastic  flow  did  occur 
such  that  W  ~  1/r  as  r  -*•  0  .  If  material  having  undergone  this  level  of 
irreversible  straining  was  left  behind  in  the  wake,  this  would  imply  that  at 
any  fixed  value  of  x^  behind  the  tip  (i.e.,  x^  <  0)  W  -  l/Jx^l  as 
Xj  0  .  But  this  is  clearly  not  possible  since  it  implies  infinite  energy 
in  any  finite  length  of  the  wake.  The  argument  is  not  rigorous  since  it  is 
conceivable,  although  it  seems  unlikely,  that,  once  plastic  strr  "ing  occurs 
such  that  W  ~  1/r  ,  this  level  of  plastic  straining  could  be  reversed  so 
that  the  residual  plastic  strains  in  the  wake  were  of  smaller  order  of 
magnitude  so  that  -*■  0  as  Xj  0  •  for  fixed  . 


where  y  is  the  surface  energy  per  unit  area  required  to  create  a  new  surface. 

The  surface  energy  y  is  introduced  into  (9)  on  physical  and  not  mathematical 
grounds.  Furthermore,  Cherepanov  goes  on  to  conclude  that,  since  Y  i  0  »  w  must 
have  a  1/r  singularity  in  all  cases.  This  is  correct  for  nonlinear  elastic 
materials  as  corroborated  in  [9,  10].  Bit  for  steady  growth  of  a  line  crack  in  an 
elastic-perfectly  plastic  material,  at  leist,  the  integral  in  (9)  is  zero;  and, 
therefore,  the  proposed  fracture  criterion  can  never  be  satisfied  wir.hin  the  context 
of  this  theory.  From  a  different  approach,  Rice  [11]  has  also  shown  hat,  for  the 
small  strain  theory  of  a  line  c.-ack  in  an  elastic-perfectly  plastic  material,  the 
surface  energy  cannot  be  accounted  for  in  the  energy  balance,  whether  the  motion  is 
steady  or  not. 

The  derivation  of  [6]  suggests  that  this  equation  holds  for  nonsteady  growth. 

As  already  mentioned,  the  Integral  in  (9)  is  not  path-independent,  in  general,  unless 
steady  conditions  prevail  or  unless  the  material  is  nonlinear  elastic.  Evidently, 
Cherepanov  has  tacitly  invoked  one  or  the  other  of  these  restrictions  in  his 
derivation. 

With  the  proper  interpretation  of  W  and  T  ,  the  integral  (8)  is  path- 
independent  for  a  full  finite  strain  theory  analogous  to  the  nonlinear  elastic 
results  of  Eshelby  [8]  and  Knowles  and  Sternberg  il2]. 

3.  IMPLICATIONS  OF  THE  PATH- INDEPENDENT  INTEGRAL  FOR  SMALL  SCALE  YIELDING 

If  stable  crack  growth  occurs,  steady  conditions  in  the  vicinity  of  the  crack 
tip  will  be  approached  once  a  crack  has  propagated  a  distance  on  the  order  of  several 
times  the  active  plastic  zone  size.  If  the  plastic  zone  size  and  the  width  of  the 
wake  of  residual  plastic  strains  are  sufficiently  small  they  will  have  a  negligible 
perturbating  effect  on  the  elastic  solution  away  from  the  tip  and  shank  of  the  crack. 
In  some  region  surrounding  the  tip,  which  is  small  compared  to  the  crack  length  and 
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other  relevant  geometric  lengths  but  large  compared  to  the  active  plastic  zone  size, 
the  field  of  the  dominant  elastic  singularity  governs.  Under  steady  conditions  the 
size  of  the  active  plastic  zone  and  the  width  of  the  wake  remain  constant  and  the 
residual  plastic  strains  outside  the  active  plastic  zone  in  the  wake  are  necessarily 
independent  of  .  We  consider  the  limiting  problem  for  steady,  small  scale 
yielding  for  a  semi-infinite  crack  with  a  semi-infinite  wake  shown  in  Fig.  2a.  A 
similar  problem  formulation  is  discussed  in  [3]  for  the  anti-plane  shear  problem. 

For  fixed  0  in  the  open  range  1 9 1  <  ir 


K  5^(0) 


13  Jui  1J 


as  r  -► 


(10) 


where  K  is  the  elastic  stress  intensity  factor  and  0 (0)  are  the  0-variations 
of  the  stress  components  associated  with  the  domLnant  singularity  of  the  elastic 
solution.  In  the  wake  (i.e.,  0  -  ±ir)  as  r  -*■  00  finite  residual  stresses  and 
strains  occur  which  are  not  known  a  priori  but  depend  on  the  deformation  history  in 
the  active  plastic  zone.  In  addition  to  (10),  it  can  also  be  asserted  that  as 

-oo  for  any  fixed  x^  outside  the  wake  (i.e.,  | X2 1  >  2h)  the  stress  field 
approaches  zero.  This  follows  from  the  fact  that  far  from  the  crack  tip  either 
half  of  the  wake  can  be  viewed  as  an  infinite  strip  of  width  2h  which  has  undergone 
Xj- independent  plastic  straining  and  which  is  attached  to  an  elastic  foundation  of 
infinite  depth.  In  such  a  problem  the  stresses  and  strains  outside  the  strip  are 
not  influenced  by  the  residual  stresses  and  strains  in  the  strip,  since  for  such 
Xj- independent  deformations  the  stiffness  of  the  strip  is  zero  compared  to  that  of 
the  underlying  elastic  foundation.  By  viewing  the  wake  far  from  the  tip  in  this 
same  way,  it  also  follows  that  ^  1  +  0  and  u^  +  0  in  the  wake  as  ^  -*  -*  . 

Let  T  be  a  circular  contour  of  radius  r  centered  at  the  tip  as  In  Fig.  2a. 
c  c 

Consider  evaluating  J8  using  this  contour  and  then  let  r  -*■  “  .  The  elastic 


-10- 


field  (10)  makes  a  finite  contribution  to  J  since  for  large  ,  W  -  l/rc  , 


for  8  <  tt  ,  and  ds  =  r  d0 

c 


Furthermore,  since  the  wake  does  not  influence 


the  elastic  field  outside  it  as  r  00  and  since  a  1/r  contribution  to  W  in 

c  c 

the  wake  portion  of  the  line  integral  (8)  makes  no  contribution  as  rc  m  ,  the 

contribution  of  (10)  to  the  line  integral  decouples  from  the  contribution  of  the 

s  2  2  2 

wake.  The  contribution  of  (10)  to  J  is  K  /E  for  plane  stress  and  K  / [  (1-v  )EJ 

for  plane  strain,  where  E  is  Young's  modulus  and  V  is  Poisson's  ratio;  i.e.,  the 

well  known  energy  release  rates  per  unit  thickness  for  the  elastic  problem.  The 

contribution  from  the  wake  is 

2h 


■2  lim  ( 

c  -  J 

A 


("-Tlul,l)dX2 


(11) 


The  second  term  in  the  integrand  vanishes  since  u.  .  ■+•  0  as  discussed  above.  Let 

lfl 

WAVE  <*<inote  the  average  residual  internal  energy  density  in  the  wake  far  from  the 


tip ,  i.e., 


2h 

W  =  —  f  W  dx 

AVE  '  2h  x.  -»  J  2 

^  A 


(12) 


Combining  K  /E  (for  plane  stress)  and  (11)  using  (12)  and  the  requirement  that 
JS  -  0  gives 


K  /E  «  4h  W 


AVE 


(13) 


For  the  problem  posed,  (13)  implies  that  the  elastic  energy  release  rate, 

2 

K  /E  ,  must  equal  the  residual  internal  energy  per  unit  length  left  behind  in  the 
wake.  Equation  (13)  should  not  be  regarded  as  a  fracture  criterion  (unless  the 
value  of  4h  WAV£  associated  with  steady  fracture  conditions  happens  to  be  known 
from  other  considerations) .  It  is  a  necessary  condition  for  steady  problems  which 
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must  hold  for  any  value  of  K  ;  It  Is  the  simple  consequence  of  the  fact  that  all 
the  elastic  energy  released  by  the  extending  crack  is  deposited  in  the  wake. 

Equation  (13)  is  a  direct  demonstration  of  the  Irwin-Orowan  modification  of  Griffith's 
original  energy  balance  arguments  for  perfectly  brittle  materials  to  include  plastic 
deformation.  The  small  strain,  elastic-plastic  theory  of  a  perfectly  sharp  line 
crack  permits  no  energy  to  be  lost  at  the  crack  tip,  as  already  discussed,  and  thus 
(13)  cannot  include  a  surface  energy  contribution.  Under  most  conditions  the  surface 
energy  per  unit  length  of  propagation  of  a  metal  is  an  extremely  small  fraction  of 
the  residual  internal  energy  per  unit  length  associated  with  plastic  deformation  so 
that  its  absence  from  (13)  is  of  little  consequence,  at  least  as  far  as  the  overall 
energy  balance  is  concerned.  It  is  not  clear,  in  general,  what  importance  the 
surface  energy  plays  in  establishing  a  fracture  criterion. 

Cherepanov  [7)  has  derived  a  relation  similar  to  (13)  using  a  Dugdale  model  of 
the  plastic  zone  with  no  consideration  of  the  wake.  However,  he  includes  the  work 
done  by  the  tractions  in  the  strip  zone  in  with  the  surface  energy  to  give  an 
effective  surface  energy  per  unit  of  surface  area  formed. 

4.  PRELIMINARY  SOLUTIONS  FOR  MODELING  WAKE  OF  GROWING  CRACK 

In  Section  5  a  model  of  steady  crack  growth  in  thin  sheets  will  be  formulated. 

A  wake  of  residual  plastic  strains  and  stresses  will  be  introduced  which  extends  from 
the  crack  tip  at  xx  -  -d  to  x1  -  -»  as  depicted  in  Fig.  2b.  (For  convenience, 
the  origin  is  now  chosen  to  coincide  with  the  leading  edge  of  the  plastic  zone.) 
Solutions  used  in  representing  this  wake  are  given  in  this  section.  Since  these 
solutions  can  be  obtained  using  complex  variable  methods  which  are  now  reasonably 
well  known,  they  will  be  given  without  derivation. 

Introduce  the  complex  variable  z  *  x^  +  ix2  and  let  <J>(z)  and  41  (z)  be  the 
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Muskhellshvili  functions  for  isotropic  plane  elasticity  theory  as  usually  ierined 
so  that 

°11  +  °22  "  W+O1) 


a22  "  °11  +  21  °12  =  +  <U ) 

Uj^  +  iu2  ■  [  (l+v)/E]  [k4> -zA'  -  ijJ] 

where  K  -  (3-v)/(l+v)  for  plane  stress  and  the  bar  denotes  complex  conjugation. 

4.1  Elastic  field  surrounding  a  circular  plastic  spot 

Let  a  circular  region  of  radius  R  and  center  z  *  x°  ix°  in  an 

o  1  2 

isotropically  elastic  infinite  sheet  undergo  a  uniform  plastic  straining  . 

-  otp 

Continuity  of  tractions  and  displacements  across  the  circular  boundary  is  maintained. 
This  is  a  plane  stress  version  of  the  Eshelby  problem  [13].  The  complex  functions 
for  the  elastic  field  outside  the  circular  plastic  spot,  lz-z0l  >  &  »  are 

<[>o(z)  -  |eAc2(z-zq)”1 

(15) 

t|»o(z)  -  EA[-(c1+c2)  (z-zq)  1  +  !c2z0(z-zo)"2  +  ic2R2(z-zo)"3] 

where 

A-’r2  •  •  c2-i<EU-E22  +  21eU>  <16) 

4.2  Elastic  field  due  to  two  circular  plastic  spots  which  are  symmetrically  placed 


on  either  side  of  a  crack 


Consider  two  circular  regions  of  equal  radius  R  .  The  center  of  the  upy  jr 

spot  is  at  z  ■  x°  +  ix°  and  the  center  of  the  lower  one  is  at  z  ;  the  restriction 
o  1  L  o 

x°  >  R  is  assumed.  Along  x2  «  0  for  <  0  the  tractions  are  zero  on  the  crack 
faces.  Let  the  upper  spot  undergo  a  uniform  plastic  straining  <eii»e22’el2^  an<* 
let  the  lower  spot  undergo  <eii»e22’”el2^  *  1316  resultin8  stress  distribution  is 

symmetric  with  respect  to  x2  ■  0  .  The  solution  for  the  elastic  field  outside  the 


J 
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two  spots  is  given  by 

<Kz)  =  4>0(z)  +  4>0(z)  +  ^(z)  +  ^(z) 
iKz)  e  i^Q(z)  +  ^(z)  +  ^(z)  +  ^(z) 

where  4>q  and  ij;  are  the  functions  in  (15).  Here  the  standard  notation 
$Q(z)  =  <j>Q(z)  ,  etc.,  is  used.  The  functions  <J>^  and  are  given  by  i|^ 

and 


(17) 


4>{(z)  *  *  |[2<J>q(z)  +  z^(z) +^(z)]  -  |  EAf(z)//z 


(18) 


where 


f(z)  - 


F 


2z  (z-z  )' 
o  o 


,-c  (z+z  )  + 
1  o 


3  2  2  2-  2- 

c.(-z  +  6zz  +  3z  z  -  3z  z  -  6zz  z  +z  z  ) 
Z  o _ o _ o _ o  o _ ot>  o 

4z  (z  -z) 
o  o 


3c.R2(-5z3  -15z2z  +5z  z2  -  z3) 
2  o _ o _ o _ ' 

16(z  -z)2z2 
o  o 


(19) 


The  square  root  is  defined  such  that  Si  ■  1  with  its  branch  line  chosen  to  lie 
along  the  negative  real  axis. 

TWo  quantities  which  will  be  of  special  interest  are  the  stress  intensity  factor 
and  the  crack  opening  displacement.  Using  the  standard  definition  for  the  stress 
intensity  factor,  one  finds 


K  -  xlifQ+  (2ttXi)1/2  a22(xr0)  -  -2v5¥  EA3?e  [f(0)] 


(20) 


where  <5?e  [z]  stands  for  the  real  part  of  z  .  For  a  pair  of  spots  whose  distances 

from  the  free  crack  surface  are  small  compared  to  the  distances  from  the  tip,  i.e., 

0,1  o  | 


X2^X1^  <K  *  *  (20)  reduces  to 


3  'xjr3/2  EA 


2i/2tF 


0 

0 

2i 

rf,  Xo  +  • 

f£P  EP  \ 
{£22,C12} 

X2 

o 

1*1* 

(21) 


The  plastic  strain  components  c^2  and  e^2  have  an  inherently  smaller  influence 
on  the  stresses  at  the  tip  than  does  when  x°/|x°|  is  small. 

With  <5(x^)  =  u2(x^,0  )  -u2(x^,0  )  for  x^  <  0  ,  the  crack  opening  displacement 
is  given  by 


<S  (xx) 


r\ 

(- 


nf1/2.*e 


[f (n)]dn 


where  n  is  a  real  integration  variable.  Let  £q  be  an  arbitrary  reference  strain. 
It  is  a  simple  matter  to  reduce  (22)  to  a  nondimensional  form  for  6/(Reo)  involving 
only  the  following  quantities:  x^lx®!  ,  x°/|x°|  ,  R/x°  ,  e^/e  ,  Ej2^e0  and 
e22^Eo  *  Curves  of  6(x^)/(Req)  calculated  numerically  using  (22)  are  shown  as 
solid  line  curves  in  Fig.  3.  In  Fig.  3a  x°/|x°|  is  taken  to  be  1/10  while  in 
Fig.  3b  it  is  3/10  ;  in  both  cases  R  is  chosen  to  be  equal  to  x°  so  that  the 
spots  just  touch  the  crack  faces  as  indicated  in  the  insert.  The  curves  identified 
by  (ei]/e0 "  *  »  ei2"G22*^  were  computed  using  those  values  for  the  parameters, 
with  similar  identifications  for  the  other  cases. 

When  x°/|x°|  is  small  one  would  expect  the  displacements  in  the  vicinity  of 

the  spot  to  be  essentially  the  same  as  those  predicted  for  a  plastic  spot  situated 

near  the  free  edge  of  a  semi-infinite  sheet  (i.e.,  where  the  tractions  vanish  along 

*2  "  0  for  all  .  By  taking  the  appropriate  limit  of  (22)  as  x°  -*■  0  one  can 

show  that  this  is  the  case.  The  dashed  line  curves  in  Fig.  3  are  the  predictions  of 

this  simpler  solution  for  the  spot  of  radius  R  centered  at  z  in  a  semi-infinite 

o 

sheet  (x2  >  0)  .  The  identification  6(x^)  ■  2u2(x^,0)  has  been  made.  Note  that 
for  x2/|x°|  ■  1/10  the  two  sets  of  predictions  are  essentially  indistinguishable 
over  the  range  of  x^  shown,  and  even  for  x2/|x°j  *  3/10  the  approximation  is 
reasonably  good  except  right  at  the  crack  tip. 
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5.  SMALL  SCALE  YIELDING  KTDEL  OF  THE  EFFECT  OF  THE  WAKE  ON  STEADY 

CRACK  GROWTH  IN  THIN  SHEETS 

As  indicated  in  Fig.  2b,  the  active  plastic  zone  will  be  modeled  by  a  Dugdale 
strip  zone  extending  from  the  crack  tip  at  x1  =  -d  to  X]^  »  0  .  The  sheet  is  slit 
along  x2  »  0  in  this  zone  and  the  tractions  are  required  to  satisfy  <J22  =  Oq  and 
°12  m  ®  on  each  face.  Identify  o  with  an  effective  yield  stress  of  the  material 
in  tension.  The  length  d  of  the  zone  will  be  adjusted  so  that  the  stresses  are 
bounded  ahead  of  the  plastic  zone. 

Small  scale  yielding  is  assumed  and  thus  far  from  the  tip  on  any  radial  line 
except  0  *=  ±it  the  dominant  singularity  of  the  elastic  solution  (10)  is  approached. 
The  stress  intensity  factor  K  is  regarded  as  prescribed.  In  terms  of  the  complex 
Muskhelishvili  functions  the  far  field  is 

♦  '  “  2i -*  K/[2(2uz)1/2]  as  |z|  +  »  (23) 

No  attempt  will  be  made  to  relate  the  details  of  the  distribution  of  the  plastic 

strains  in  the  wake  to  the  deformation  field  in  the  active  plastic  zone,  but  it  will 

oe  assumed  that  all  the  elastic  energy  released  goes  into  the  internal  energy  of  the 

wake  consistent  with  (13).  In  other  words,  the  work  done  by  the  tractions  ■  a 

22  o 

in  the  strip  zone  is  considered  to  appear  as  part  of  the  internal  energy  left  in 
the  wake. 

£aS  rePresent  the  inplane  plastic  strains  averaged  with  respect  to  x2 
across  the  wake;  these  are  independent  of  x^  .  The  wake  is  modeled  by  integrating 
contributions  of  the  plastic  spot  solution  (17)  from  x£  -  -»  to  x°  -  -d  .  Identify 
the  uniform  plastic  strains  in  an  elemental  spot  with  the  average  plastic  strains  in 
the  wake  .  At  a  typical  point  x°  in  the  wake  identify  the  area  element  A 

in  (16)  and  (18)  with  2h  dx?  and  take  the  element  to  be  centered  at  z  *  x?  +  ih 
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with  its  partner  centered  at  zq  .  Note  that  (17)  makes  no  contribution  to  the 
tractions  on  the  surfaces  in  the  strip  zone  ahead  of  the  crack. 

The  energy  relationship,  K2/E  =  4h  WAV£  ,  discussed  in  §3  will  be  imposed.  It 
is  convenient  to  introduce  a  set  of  dimensionless  quantities  g  D  defined  by  the 

Ctp  J 

equation 


aoea6  =  gaB  WAVE 


These  quantities  reflect  the  relative  amounts  of  each  component  of  the  average  strain 

in  the  wake.  For  example,  if  is  the  only  nonzero  inplane  component  generated 

in  the  active  zone,  then  the  average  residual  internal  energy  in  the  wake,  WA1M_  ,  is 

^  AVE 

equal  to  ^  material  is  elastic-perfectly  plastic  and,  thus,  g  ^  m  1  and 

8ll“812“°  *  In  general,  the  elastic  residual  energy  in  the  wake  will  comprise  only 
a  relatively  small  fraction  of  WAy£  so  that  the  quantities  g^  will  range  in 
magnitude  from  zero  to  approximately  unity. 

At  this  stage  the  width  of  the  wake,  4h  ,  will  be  left  unspecified  except  that 
it  will  be  assumed  that  the  ratio  h/d  is  small  compared  to  unity,  not  exceeding 
about  4/10  ,  consistent  with  theoretical  and  experimental  results  on  the  shape  of 
plane  stress  plastic  zones.  Later  it  will  be  seen  that  the  main  predictions  of  the 
model  are  essentially  independent  of  the  choice  of  h  . 


The  amplitude  of  the  inverse  square  root  singularity  of  C^Cx^O)  ahead  of 
the  active  plastic  zone  is  the  sum  of  the  contributions  from  (23),  from  the  traction 
°22  m  °o  applied  along  the  strip  zone  and  from  the  wake.  The  first  two  are  the 
contributions  which  enter  into  the  Dugdale-Barenblatt  model  and  the  third  is  obtained 
by  integrating  the  contributions  in  (20)  over  the  entire  wake.  Thus, 


'2ttxi  a22(x1»0)  -  K  -  Oq 


-a 

(8d/7T) 1/2  -  4>/2tt  Eh  J  dx°  iSfe  (f  (0)  ] 


J 


njHTillHtW 


The  radius  R  of  the  spot  appears  explicitly  only  in  the  third  term  in  the 
brackets  in  the  expression  (19)  for  f(z)  .  For  the  time  being  we  will  retain  the 
term  in  (19)  which  explicitly  involves  R  without  identifying  it  with  a  specific 
length,  but  it  will  be  assumed  that  R  is  not  greater  than  h  .  Later  it  will  be 
seen  that  this  term  has  very  small  influence  on  the  predictions. 

For  the  stationary  crack  the  third  term  in  (25)  is  absent  and  d  is  chosen 
such  that  there  is  no  inverse  square  root  singularity  in  the  stress  O ^  at  the 
leading  edge  of  the  plastic  zone;  i.e.,  for  the  Dugdale-Barenblatt  model. 


dD  *  (u/8)(K/ao)2  (26) 

The  third  term  in  (25)  can  be  integrated  with  the  result  that  this  term  is  (letting 
b  =  -d  +  ih) 

-4^2?  Eh.^e  [  (Cj+c^b-172  -  -  ih  c2b‘3/2  -  ~  c2rV5/2] 

Expanding  this  expression  in  a  Taylor  series  in  small  values  of  h/d  gives 

Eh  d“1/2iJ?e  [  (h/d)  (Cj+2c2)  +  0{h/d)2] 

-  -E(2d/7i)1/2(h/d)2[eJ1  +  £>{(h/d)e^2>  (h/d)e22}] 


Now,  use  (13),  (24)  and  (26)  to  rewrite  the  above  expression  for  the  contribution 
of  the  wake  to  (25)  in  terms  of  the  8ag's  with  the  result 


(27) 


For  the  steady  growing  crack  d  will  also  be  chosen  such  that  (25)  vanishes. 
The  underlined  term  in  (27)  will  be  only  a  few  percent  if  h/d  is  not  larger  than 
about  3/10  and  if  is  not  the  predominant  plastic  strain  component  so  that 

the  magnitude  of  g^  will  not  exceed,  say,  0.3  to  0.4  .  Thus,  the  contribution 
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of  the  wake  (27)  will  not  exceed  a  few  percent  of  the  second  term  in  (25).  Therefore, 
the  length  of  the  active  plastic  zone  for  the  steadily  growing  crack  is  almost 
identical  to  that  of  the  stationary  crack  at  the  same  value  of  K  ,  i.e.,  d  =  dp  . 

This  result  is  consistent  with  the  similar  finding  of  Chitaley  and  McClintock  [3]  in 
their  complete  analysis  of  the  anti-plane  shear  problem. 

The  total  opening  displacement  in  the  active  plastic  zone  is  the  sum  of  the  same 
three  contributions  mentioned  above.  Using  (22)  for  the  contribution  from  the  wake 
plus  the  two  contributions  present  in  the  stationary  strip  model  gives,  for 
-d  <  xx  <  0  , 


<5^) 


4_  _o 
TT  E 


X1 
1+  d 


£n 


1+  (-Xl/d) 


1/2 


1-  (-x^/ d) 


1/2 


+  16h 


—OO 


|  dn(-n)-1/2,>?e  [f(n)l 

o 


(28) 


To  calculate  6(x^)  for  a  prescribed  value  of  K  it  is  necessary  to  first 
choose  d  such  that  (25)  vanishes.  We  will  simplify  this  procedure  by  taking 
d  *  dp  ,  which  as  already  discussed  is  an  excellent  approximation  if  h/d  is  not 
too  large  and  is  exact  for  h/d  ■+■  0  .  The  integrations  in  (28)  cannot  be  carried 
out  in  closed  form,  but  an  excellent  approximation  can  be  obtained  for  sm?»ll  h/d 
using  the  simpler  solution  alluded  to  in  §4.2  for  the  spot  in  a  semi- inf inite  sheet. 
The  approximate  expression  for  6(x^)  is  (with  d  *  dp  and  for  -d  <  x^  <  0) 


«(x1) 


8  V 

IT  E 


1/2 


1  +  • 


'  rl+(-x1/d)1/2'1 

jin'  1 


>1/2 


l_l  -  (-Xj^/d) 

+  16h  ,7m  j(c1+c2)Jtn(l  -x^b)  +  ih  c2 j^b-Xj)"1  -  b-1j  +  j  C2R2 j^b-x^2  -  b  2  | 


(29) 
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where  .9 m  [z]  denotes  the  imaginary  part  of  z  ,  b  =  -d  +  ih  ,  and  £n[z]  is 
the  natural  logarithm  of  complex  argument  defined  such  that  Jin [ 1 ]  =  0  with  the 
branch  line  taken  along  the  negative  real  axis. 

Let  <5t  =  6(x^  =  -d)  denote  the  crack  tip  opening  displacement.  By  setting 
x^  =»  -d  in  (29)  and  rearranging  the  result,  one  can  obtain 

<$t  -  (8/TT)Ood/E  -  8TTh(C;L+c2)  +  4h[4-  (R/h)2]^m  [c2] 

+  8h  |-2(Cl+c2)in(-b/d)  -  2ih  c^"1  -  ~(R/h)2c2h2b~2  j-  (30) 

Denote  by  6^  the  crack  tip  opening  displacement  for  the  stationary  crack;  i.e., 
the  Dugdale-Barenblatt  value 

6®  *  (8/7 r)a  d/E  =  K 2/(a  E)  (31) 

to  o 

Here  again,  use  the  relation  (13)  between  K  and  the  residual  internal  energy  in 
the  wake  together  with  (24),  and  (26)  to  express  5t  in  terms  of  the  g^'s  with 
the  result 

V6t  "  1  '  l(8ll+822)  +  (R/h)2^i2  +  ^(h/d)  (32) 

The  terms  indicated  by  £?(h/d)  are  the  terms  from  the  second  line  of  (30);  they 
go  to  zero  linearly  with  h/d  as  h/d  -*■  0  .  For  the  mathematical  limit  h/d  0  , 
(32)  gives  the  exact  expression  for  as  based  on  the  unapproximated  equations 

(25)  and  (28). 

For  small  h/d  ,  as  given  by  (32)  is  essentially  independent  of  the  choice 
of  the  width  4h  of  the  wake  by  virtue  of  the  fact  that  the  quantities  ga^ 
introduced  in  (24)  are  themselves  independent  of  h  ,  except  to  the  extent  that 
the  relative  proportions  of  the  average  plastic  strain  components  in  the  wake  may 
be  a  function  of  h/d  .  Note  that  the  term  involving  R  ,  which  has  been  retained 
in  the  analysis,  enters  only  in  the  one  term.  Its  influence  is  clearly  small  since 


J 
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over  its  meaningful  range,  0  <  R  <  h  ,  the  coefficient  of  g^  in  (32)  changes  by 
only  twenty-five  percent. 


The  values  of  g^^  depend  on  the  details  of  the  plastic  deformation 
occurring  in  the  active  plastic  zone.  For  the  purpose  of  discussion  suppose 
that  e^2  is  the  predominant  inplane  plastic  strain  component  so  that  g00  =  1 


22 


and  =  g^2  "  0  .  In  this  case,  according  to  (32),  the  crack  tip  opening 

displacement  of  the  steadily  growing  crack  is  about  one-half  the  value  of  the 
stationary  crack  at  the  same  value  of  K  .  If  5  is  used  as  a  measure  of  the 
intensity  of  deformation  at  the  crack  tip  and  if  a  fracture  criterion  is  used 
which  involves  the  attainment  of  a  critical  value  of  6  ,  then  the  model 

indicates  that,  due  to  the  wake  alone,  the  critical  value  of  K  for  steady  growth 
should  be  approximately  times  the  value  for  initiation. 


The  fact  that  the  wake  reduces  the  crack  tip  displacement  by  an  amount  6^/2 
when  £^2  is  predominant  can  be  seen  very  simply  in  the  following  way.  Neglect 
the  interaction  with  the  end  of  the  zone  at  x^  ■  0  and  consider  a  semi- infinite 
half  sheet  with  a  traction-free  edge  on  X2  *  0  .  Let  a  uniform  plastic  straining 
e^2  be  induced  in  the  wake  region,  0  <  x2  <  2h  and  x^  <  -d  .  The  displacement 
u„  for  this  problem  .at  .  (x^*  -d,x2 ■  0)  is  exactly  -he£„  .  Using  V22  *  WAVE 
and  (13)  to  express  -he?2  In  terms  of  K  gives 

u2(-d,0)  -  -K2/(4aQE)  (33) 

Since  the  total  contribution  to  6  from  both  halves  of  the  wake  is  twice  the  above 
amount,  it  is  immediately  seen  from  (31)  that  this  is  exactly  -6^/2  .  This  same 
result  (33)  holds  for  any  distribution  of  e22^X2^  in  the  walce  reSion  as  lon8  as 
it  is  independent  of  x^  .  This  observation  lends  some  confidence  in  our 
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representation  of  the  plastic  strains  in  the  wake  as  being  essentially  uniform. 

Generally  and  will  also  be  present  in  the  wake  and  their 

contributions  can  also  be  expected  to  diminish  6^  relative  to  6^  since  various 

t  t 

evidence  points  to  positive  and  negative  g^  •  In  particular,  the  trailing 

edges  of  an  advancing  crack  in  a  thin  sheet  are  often  observed  to  buckle  out  of  the 
plane  which  is  most  likely  due  to  a  residual  compressive  stress  0^  associated 
with  positive  .  If  lines  are  lightly  scored  perpendicular  to  the  line  of  the 

advancing  crack  in  the  face  of  the  sheet  ahead  of  the  tip,  following  passage  of 
the  crack  the  lines  indicate  negative  values  of  e^2  above  the  crack  and  positive 
values  below,  consistent  with  negative  g ^  •  The  value  6 1 / 6^  =  1/2  seems  to  be 
more-or-less  representative  for  any  such  combination  of  g^’s  constrained  by  (24). 

The  result  (32)  holds  for  plane  strain  as  well  as  plane  stress,  although  the 
strip  model  of  the  plastic  zone  is  not  realistic  in  plane  strain.  Nevertheless, 

(32)  is  very  suggestive  as  far  as  the  effect  of  the  wake  is  concerned.  In  plane 
strain  is  small  and,  since  the  plastic  volume  change  Is  zero,  this  implies 

that  +  e^2  s  0  or  g^  +  g22  =  0  .  Thus,  the  term  g22) /2  in  (32), 

which  gives  the  major  contribution  to  reducing  the  crack  tip  opening  displacement 
in  plane  stress,  is  negligible  in  plane  strain.  This  suggests  that  the  effect  of 
the  wake  on  stabilizing  crack  growth  is  less  in  plane  strain  than  in  plane  stress. 

Figure  4  shows  curves  of  5(x^)  in  the  strip  zone  for  three  combinations  of 
g^g's  from  which  any  other  combination  can  be  computed  at  the  same  value  of  h/d 
by  linear  superposition.  The  curves  were  calculated  using  (29)  for  h/d  =  0.2 
and  R  ■  h  .  The  smaller  is  h/d  ,  the  more  the  effect  of  the  wake  on  6(x^)  in 
the  strip  zone  becomes  concentrated  near  x^  *  -d  .  Thus,  for  example,  for  g22  *  1 
and  g,,  ■  g,  „  *  0  ,  6^  is  no  longer  the  maximum  value  of  S(x-)  on  the  interval 
-d  <  x^  <  0  when  h/d  is  less  than  about  0.2  .  The  model  as  formulated  can  no 
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longer  be  considered  physically  relevant  when  S(x^)  does  not  decrease  monotonicalJ y 

from  the  crack  tip  to  the  leading  edge  of  the  strip  zone.  The  cutoff  value  of  h/d 

ranges  from  approximately  .1  to  .2  depending  on  the  combination  of  g  's.  Thus. 

ap 

while  (32)  is  exact  for  the  mathematical  limit  h/d  ->  0  ,  the  model  loses  its 
significance  before  this  limit  is  reached.  In  particular,  its  relevance  to  stable 
propagation  in  very  thin  sheets  becomes  questionable  when  the  plastic  deformation  is 
largely  confined  to  a  thin  zone  of  width  comparable  to  the  sheet  thickness  in  which 
necking  occurs.  Then,  if  the  sheet  is  very  thin  the  length  of  the  active  zone  may 
be  many  times  the  sheet  thickness  such  that  h/d  falls  outside  the  range  of  validity 
of  the  model. 

We  reiterate  that  irreversibility  associated  with  nonproportional  loading  in 
the  active  plastic  zone,  which  also  contributes  to  a  stabilization  of  crack  growth, 
has  not  been  taken  into  account  in  this  study.  A  complete  analysis  would  necessarily 
involve  an  accurate  modeling  of  both  the  active  plastic  zone  and  the  wake  such  as 
that  carried  out  in  [3]  (but,  perhaps,  based  on  a  stress-strain  relation  with  a 
non-smooth  yield  surface).  The  main  conclusion  of  this  study  is  that  the  wake  of 
residual  stresses  and  strains  contributes  a  significant  stabilizing  effect  to  crack 
growth. 
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FIG.  3  CRACK  OPENING  DISPLACEMENT  RESULTING  FROM  A  PAIR 
OF  CIRCULAR  SPOTS  WHICH  HAVE  UNDERGONE  UNIFORM 
PLASTIC  STRAINING 


STRIP  ZONE 


OPENING  DISPLACEMENT  IN  STRIP  ZONE  FOR  STEADILY 
GROWING  CRACK  COMPARED  WITH  THAT  OF  STATIONARY 
CRACK  (h/d  =  2/10). 


